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1.  Introduction 


Groundwave  propagation  is  one  of  the  oldest  studied  topics  in  ap¬ 
plied  electromagnetics.  Several  models  of  the  propagation  medium  have 
been  considered,  e.g.,  flat  or  spherical  Earth,  homogeneous,  stratified 
or  anisotropic  ground,  corrugated  surfaces,  single  or  mixed  path  propa¬ 
gation.  An  excellent  summary  of  problems  and  modern  solutions  to  this 
subject  is  given  in  [lj  ,  which  contains  also  a  large  number  of  referen¬ 
ces  . 

A  canonical  problem  in  mixed-path  propagation  is  that  of  an  abrupt 
discontinuity  in  ground  parameters  along  a  straight  line,  as  depicted 
in  Fig.  1;  the  ground  exibits  different  properties,  e.g.,  conductivity 
and/or  permittivity,  for  y  >  0  and  y  <  0.  The  canonical  problem  is  the 
following:  for  a  given  incident  plane  wave,  compute  the  field  at  point 
P.  If  the  plane  wave  solution  is  needed  to  synthesize  the  field  produced 
by  localized  sources,  the  angle  <l>0  shpuld  be  allowed  to  be  complex,  and 
the  incidence  should  not  be  restricted  to  the  normal  case  ( two-dimensio_ 
nal  problem) . 

A  central  role  in  the  solution  to  this  problem  is  played  by  the 
surface  impedance  concept,  i.e.,  by  definition,  the  ratio  of  the  tangen¬ 
tial  components,  and  H  of  the  fields  at  the  interface  x  -  0.  For  the 

O  ti 
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problem  of  an  incident  plane  wave  with  a  fixed  angle  <p0,  the  use  of  an 
impedance  type  boundary  condition  is  an  exact  approach.  However,  when 
using  plane  wave  expansions  (as  in  the  solution  of  the  just  mentioned 
canonical  problem) ,  the  use  of  an  angle  independent  surface  impedance 
implies  an  approximation  to  the  exact  boundary  conditions. 

A  discussion  about  the  validity  of  the  surface  impedance  concept 
is  given  in  [2]  and,  more  recently,  in  [3]  .  We  want  also  to  mention 
a  few  of  the  excellent  agreements  between  experimental  and  theoretical 
results  using  impedance  type  boundary  conditions:  see  [U,  5]  for  the 
case  of  scattering  by  a  half-plane  with  two  face  impedances,  and [6-8] 
for  the  case  of  laboratory  models  of  groundwave  propagation.  Additional 
results  are  listed  in  [l]  .  Accordingly,  the  use  of  impedance  type  boun¬ 
dary  conditions  seems  quite  adequate  to  the  case  at  hand. 

The  solution  to  the  canonical  problem  depicted  in  Fig.  1  is  known 
essentially  for  the  case  of  normal  incidence  [9~l6]  ;  when  oblique  inci¬ 
dence  is  considered,  only  an  approximate  matching  of  tangential  compo¬ 
nents  of  field  at  the  impedance  discontinuity  line  is  provided;  and 
restrictions  on  surface  impedances  are  imposed,  e.g.,  one  of  them  is  ze¬ 
ro  [lT-18]  .  We  note  that  the  problem  of  scattering,  under  normal  inci- 
ence,  by  a  wedge  of  arbitrary  angle  with  two  different  face  impedances 
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has  been  rigorously  solved  since  1958  by  Maliuzhinets  [  19  ]  .  When  the 
wedge  angle  is  made  equal  to  it,  the  problems  of  Fig.  1  is  recovered. 
Maliuzhinets '  solution  has  been  recently  extended  to  the  case  of  obli¬ 
que  incidence,  for  the  half-plane  case  [ 20]  ,  and  to  the  case  of  a  pla¬ 
nar  surface  impedance  discontinuity  [21]  . 

Before  considering  the  extension  of  Maliuzhinets’  solution  and 
its  modifications  appropriate  to  groundwave  propagation  problems,  it 
is  worth  remembering  that  an  alternative  convenient  way  to  the  study 
mixed-path  propagation  problems  is  obtained  by  using  the  compensation 
theorem  [22]  .  This  procedure  has  been  originally  exploited  by  Wait 
[23-24]  and  then  applied  to  oblique  propagation  of  groundwaves  across 
a  coastline  [25-27]  .  Extensive  theoretical  and  experimental  work, 
using  essentially  this  procedure,  is  referred  to  in[  1  ] .  Use  of  the 
compensation  theorem  is  very  attractive  inasmuch  as  it  allows  one  to 
consider  not  only  the  simple  canonical  problem  of  Fig.  1,  but,  in  prin¬ 
ciple,  any  geometry  of  surface  discontinuity,  therefore  rendering  the 
procedure  very  powerful  from  the  application  viewpoint.  However,  the 
resulting  integral  equation  should  be  solved  numerically  in  general, 
even  when  some  simplifying  approximations  are  made.  Accordingly  the  ri¬ 


gorous  solution  of  the  canonical  problem,  e.g.,  the  one  depicted  in 


-  8  - 


Fig.  i,  is  worth  inasmch  a  physical  description  of  the  refraction  pheno¬ 
mena  at  the  impedance  discontinuity  is  obtained.  We  also  mention  that 
mode-matching  techniques  have  been  used  for  solving  idealized  canonical 
problems  [28]  .  Comparison  between  different  theories  is  given  in  Ref. 
|29|  . 

In  the  rigorous  solution  of  the  problem  of  Fig.  1  a  key  role  is 

+  — 

obviously  played  by  the  two  surface  impedances  Z  and  Z  or,  better,  by 
the  related  parameters 


■  n±  Z 

sm  0  = 

e;  sin  0O 


(1.  1) 


where  B0  is  related  to  the  ablique  incidence  angle  (see  Fig.  2  and 
1.  1)  and  the  ±  signs  refer  to  y  >  0  (plus)  and  y  <  0  (minus)  respecti¬ 
vely  (when  confusion  does  not  arise,  we  omit  the  superscripts). 

For  groundwave  propagation 


wherein  e  is  the  relative  permeability,  a  the  conductivity  and  6  the 
r 

skin  dept^;  a  time  dependance  exp(  -iwt)  has  been  assumed;  and  the  appro¬ 
ximation  requires: 


u>e0e  «  a 

°  r 


(1.  3) 
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Inequality  (1.  3)  is  satisfied  up  to  frequencies  f  -  10kHz  for 
-1+ 

a  very  dry  soil  (a  =  10  siemens /m;  z  -  15)  and  up  f  -  1  MHz  for 

r 

-2 

normal  soil  (a  -  10  siemens/m,  z  =15);  over  the  sea  (a  -  5  sie- 

r 

mens/m,  z  -  80)  the  inequality  is  valid  in  all  the  frequency  range 
r 

of  groundwave  propagation.  In  any  case,  the  (complex)  angle  ©defined 
by  (1.  1)  is  always  small,  and  usually  very  small.  We  are  consisten¬ 
tly  using  this  last  assumption  in  all  the  body  of  this  report.  Accor¬ 
dingly,  we  take  (for  0  =  90°) 

e*  |  »  “(1-i),  [e]  «  1  (1.  U) 

It  is  also  necessary  to  introduce  another  characteristic  angle,  rela¬ 
ted  to  the  normalized  surface  admittance,  hence: 

sin  y~  =  — ~ ^  (1.  5) 

Z-  sing0 

The  (complex)  angle  y  can  be  easily  computed  under  the  approximation 
(1.  2),  hence  (for  0  =  90°) 


y  *  l + 1  ln  7r 3  u " 1  ln 

[y]>  1 


(1.  6) 
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In  passing,  we  note  that  the  underlying  assumptions  leading  to 
(1.  4  and  6)  are  also  sufficient  to  validate  the  use  of  the  surface  impe¬ 
dance  concept  for  non  plane  wave  excitation. 

It  has  been  first  recognized  by  Senior  [l7~l8]  that  propagation  of 
waves  along  a  mixed-path,  as  depicted  in  Fig.  1,  consist  assentially  of  a 
scattering  problem,  although  it  is  usually  referred  to  as  a  refraction 
one,  in  the  frame  of  surface  propagation  (for  istance,  costal  refraction). 
This  is  clearly  shown  in  the  Maluizhinets '  theory,  wherein  the  total 
field  is  decomposed  in  its  various  components:  the  incident  one  (a  pla¬ 
ne  wave),  the  two  reflected  waves  on  the  two  surfaces  y  <  0  and  y  >  0, 
the  two  surface  waves  along  y  <  0  and  y  >  0,  and  the  (cylindrical )  wave 
scattered  by  the  line  impedance  discontinuity  y  =  0.  Particularly  impor¬ 
tant  are  the  surface  waves,  since  they  are  characteristic  of  the  surface 
and,  therefore,  less  sensitive  to  the  approximation  of  the  model,  e.g., 
undulation  of  the  surface  or  Earth's  curvature. 

In  the  two-dimensional  case  of  Fig.  1  ( z-independent  fields)  an 
H-polarized  surface  field  (the  only  of  interest  in  radiopropagation  along 
the  Earth's  surface)  is  characterized  by  a  propagation  factor (the  wave 
progresses  along  -  y); 

exp(-ikx  sin  6-iky  cos  6  )=  exp  [-ikp  cos(<J>-  ~  +6  )]  (1.  7) 

2  2. 

wherein  k  =  a)  e0y0  and  is  the  free-space  propagation  constant.  According 


-li¬ 


the  incident  angle  is 


(1.  8) 


For  the  three-dimensional  case  of  Fig.  3,  wherein  the  field  propaga¬ 
tes  with  a  (real)  angle  8  with  respect  to  the  z-axis,  the  appropriate 
propagation  factor  is  obtained  from  (l.  7)  by  using  the  rotation  of  the 
coordinates 

y  =  y  sin  8  -  z  cos  8  (1.  9) 

Accordingly,  we  get  the  new  propagation  factor 


exp(-ikx  sin  9 -iky  sin8  cos  0+  ikz  cos8  cos  0  )= 

(1.  10) 

=  exp[~ik  p  cos (41  -  <f>  0)]  exp(ikz  cos80) 

Therefore,  the  incidence  angles  <f>0,  80  are  given  by 


cos 

80 

=  cos 

0  cos8 

sin  9 

cos 

*0 

’/I 

2  2 
-  cos  8  cos  8 

cos  8  sinB 

sin 

$0 

'  f 

2  2  _ 

C.  C.Q 

COS  0  COS  P 


(1.  11) 


Accordingly,  to  the  first  order  in  9 


(1.  12) 


The  simple  source  model  for  exciting  surface  waves  over  an  impedance 
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plane  is  an  infinite  magnetic  line  of  intensity  I  ,  as  depicted  in 

m 

Fig.  i.The  free  space  magnetic  field  produced  by  the  line  is  z-polari- 
zed  and  given  by 
k  I 


H  =  - —  H0(l)(kd) 

k  C 


(1.  13) 


wherein  (x)  is  the  Hankel  function  of  the  first  kind.  Let  us  use 

the  integral  representation  of  the  Hankel  function: 


,(D 


Hr'Ud)  =  ± 


i  r 


exp  [ikd  cos(a-n  )]  d  a  = 


1 

IT 


's  (1.  11+) 

exp[ikPo  cos(a-  ri0)]  exp[  -ikP  cos($-  "T  +  a)]  da 


wherein  S  in  the  Sommerfield  contour  of  integration  depicted  in  Fig. 
1+.  Inspection  of  (1.  ll+  and  13)  shows  that  the  line  source  is  synthe¬ 
sized  as  superposition  of  plane  waves  of  (variable)  incidence  angle 


°  ~  ~2  ~  a 


(1.  15) 


and  of  spectral  intensity 


exp  fikp  cos(a-  f)0)] 


(1.  16) 


Accordingly,  superposition  of  the  plane  wave  solution  weighted  by  the 
factor  (1.  l6)  allows  the  computation  of  the  field  produced  by  the 
line  source . 


As  a  last  remark  in  this  introductory  section,  let  us  state  that  the 
Sommerfeld  contour  of  integration  S,  after  deformation  into  the  steep¬ 
est  descent  path,  transforms  in  the  real  axis,  from  -  “  up  to  +  “  , 
upon  use  of  the  variable  change 

/F exp( i  ^  )  sin  =  t  (l,  17) 


2.  Mathematical  background 

In  the  Muliuzhinets *  theory  of  scattering .by  the  configuration 
depicted  in  Fig.  1,  a  central  role  in  played  by  the  Maliuzhinets '  fun¬ 
ctions  [19] 

V.  «?(a,  6)  ;T*  T  (a,y  )  (2.1) 

h  h  e  e 

The  two  functions  are  identical,  the  second  being  given  by  the  first 
one  upon  replacement  of  Y  with©  .  Accordingly,  only  the  function  V 
is  presented  hereafter,  and  the  subscript  and  0  dependance  is  dropped. 
We  have: 

¥(a)  =  N(a  +TT  -  G)  N  (a  +  Q  )  N  (a  -tt  +  ©”)  N  (a  -  ©  )  (2.  2) 


Or,  equivalently: 

*(«)  -  N^(~)  {I 

2  N(a  -  Q)  N  (a  +  ©  ) 


«  -e  ■»  - 


cos 


a  -f  — 

2  , 

cos  “  (2. 


The  function  N  (v)  is  given  by 


N  ( v ) 


2t  -  Tisint 
cos  t 


dt 


(2.  k  ) 


and  is  a  oarticular  determination  of  a  most  general  class  of 
function  introduced  by  Maliuzhinets  [19]  with  reference  to 

the  problem  of  impedance  wedge  scattering.  The  function 

N(v)  is  easy  to  compute  numerically  and,  in  any  case,  it  is 

simply  related  to  a  tabulated  function  [20,  30]  .  Its  main 

properties  are  hereafter  summarized: 


N ( V )  =  N ( -V  ) 


N  ( V  +  -)  N  ( v  -  -)  *  N2(~)  cos  - 
2  2  2  2 


lim  N( v) 


r  Lvji 

[exP  —  J  5 


v  =  v ' +  i  v" 


(2.  5) 

(2.  6) 

(2.  7) 


The  location  of  poles  of  'H(a)  within  the  strip 


3tt 


<  a'  < 


3it 


is  of  interest.  From  (2.  4)  it  is  noted  that 
a 

2t  ir  sin  t 
cos  t 


dt 


(2.  8) 


is  finite  for  a=7T/2,  diverges  toward  negative  values  for 
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a  *3tt/2  and  toward  positive  values  for  a  =  5  tt  /  2 .  Since  the 
divergence  of  (2.  8)  is  of  logarithmic  type,  the  function 
N(a)has  a  simple  pole  at<a=  ±5^/2.  Accordingly,  the  pole  of 
¥('a)  in  the  strip  of  interest  are,  from  (2.  2  ): 


«*  •  ±<  f  *  eJ) 


(2.  9 ) 


An  appropriate  use  of  relations  (2.2  -  3  -  5  -  6)  may 
avoid  the  necessity  of  effective  computation  of  the  function 
^(a).  We  quote  here  after  those  cases  which  are  relevant  to 
subsequent  analysis. 

Use  of  (2.  3  -  5  ~  6)  gives 

N2(f)N(J  +e+  +  t)H(J  -  0_  -  T  )K(f  -  -e'-T) 

Va +  -  — 1 — - 1 1 - 


n(  \  -  e++  T  ) 


e  -x 
•  sm  — - - 


(2.  10) 


If,  as  assumed,  0  is  small  and  t  is  prescribed  to  assume  on¬ 
ly  small  values ,  then 


...  ,  IT  \  „ 4 ,  7T  .  0+-T 

V  5  *  T)  ‘  *  (  i  >  ,ln  ~i~ 


(2.  11) 


Similarly 


3tt  k  xr  N(2  "  T  +  9+)N(  2  ~  T  "  0_) 

v*  V +-  t)  ■  ®  {v - 

n (  ~-t-&  +  )n(  2~t  +  6  ^ 


1 


a  i  ‘  mm 
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8  +  T  6  -T 

sin  - “ -  sin  - - - 

2  2 


e  +  x 


..U,^  r  .  e  +x  .  e  -t  /  .  e 

:N  (-)  sm  - - -  sin  - “  /sin  — 


under  the  same  previous  assumptions . 


+  T 

> 


(2.  12) 


To  complete  this  section  let  us  q.uote  an  integral  which 
is  systematically  used  in  the  body  of  this  Report  and  which 
can  be  evaluated  properly  accomodating  results  of  Ref.[3l]: 


I(ft)  =  /2  exp  ( i-j^  +  i  ft  )  sin  ^ 


8  exp(-  flt‘ 


2  .  .2  0 

J  t  -  2i  sin  ~~~ 

-CO  2 


dt  = 


F  {  /ifF  sin  ~) 


=  iri  exp(  i  fl  cos  0  ) 


F  (0) 


(2.  13) 


wherein  F(x)  is  the  Fresnel  integral 


F(x)  *  exp(i  t2)dt 


(2.  lU) 


3.  Mixed-path  -propagation .  Formal  solution  for  plane  wave 


oblique  incidence 


Let  us  consider  the  (three-dimensional)  problem  of  plane 
wave  scattering  by  the  mixed-path  configuration  depicted  in 
Fig.  2.  Any  component  F  the  incident  plane  wave  is  represented 
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by 

G1(P»  4>.z)  =  Go  exp(ik.2  cos  Q  0  )  exp  [-ik  p  sin£0  cos  ( <p~<t>0  )  ] 

(3.  1) 

Without  loss  of  generality  we  take  0  <  Re  [j30]_<  ~  and 

drop  the  z-dependence  of  the  field  components • 

Boundary  conditions  are  expressed  as 

E-i  i  •  E  =  Z”  i  xH  a  t  <f)  =  ±7  (3.2) 

—  —  x  —X  ~  -X  —  2 

or,  equivalently, 

H-ii  H  =  -7.  E  x  i  at  <f>®  ±~  (3.3) 

~  ~x~x  ~  2-  ~  —  x  2 

By  taking  the  divergence  of  (3.  2  and  3)  and  using  Maxwell's 

equations,  we  get 


SE  + 

—  7—7  +  ik  sin  £"  E 

p  6  <j)  x 


at  <j>=±‘ 


—  ~ — 7  +  ik  sin  y  H 

p  <S  <f  x 


,  IT 

at  <j>=±~ 


(3.  k) 

(3.  5) 


respectively . 

Let  us  now  expand  any  component  G  of  the  total  field  as 

a  superposition  of  plane  waves,  hence: 

’  * 

G(p,4>)  =  ~T  G(a+<|>  )exp(-ikP  s  in(30cosa )  da  (3.  6) 

}T 

where  T  is  the  two-branch  contour  of  integration  depicted  in 
Fig.  5. 

Forcing  boundary  conditions  (3.  **  and  5)  in  (3.  6)  we 


get,  upon  integration  by  parts 


(sina±sinG  )E  (  a±~)exp{ -ik  p  si  nB  0  cosa  )  da=  0  (  3  •  7) 
j  x  2 

r 

t 

—  ^  TT 

(sinaisiny  )H  ( a±— )exp( -ikp  s i n0ocosa  )  da^  0  ( 3 .  8) 
)  x  2 

r 


Solution  to  equations  of  type  {3.  1  ~  8)  is  given  inj[l9] 


under  the  assumtion  that  a  simple  pole  does  exit  at  = 


G(o)  * 


Gc  cos  $0 
sin  a-  sin 


+  E  G  s  i  n  a 
„n  n 


JLIslL- 

¥U0) 


(3.  9) 


wherein  G^  is  clearly  the  residue  at  <J>  =  <j>0  ,  G^  are  arbitrary 
constants  and  the  function  fCct)  is  the  Maliuznets'  function 
introduced  in  Sect.  2. 

In  order  to  determine  these  unknown  constants  G  let  us 

n 

express  the  spectral  z-components  of  the  field  (E  ,  H  )  in 

z  z 

terms  of  the  x-components  (E  ,  H  )  .  We  have  (see  Appendix  A) 

x  x 

cosacosBE  (a)-sina£H  (a) 

E  (a)  =  - - - - - sinB0  (3.  10) 

2  2  2 

1-cos  asin  B0 

a  ^ 

cosa  cosBCH  (a)-sinaE  (a) 

CH  (a)  *• - 2 - 2 —  sin  B0  (3-  11) 

~ Z  2  2 

1-cos  asin  B0 


As  discussed  in  [32],  the  value  of  any  field  components 
G  for  p-*-0  is  related  to  the  asymptotic  behaviour  of  the  corre¬ 
sponding  spectral  component  G  as  [a"]  -*■<*>.  In  particular,  G  should 
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be  bounded  if  G  is  required  to  be  regular  as  it  happens  for 

E  and  H  . 
z  z 

Substituting  the  expansion  (3.  9)  in  (3.  10  and  ll) 

and  using  (2.  7),  it  is  easy  to  convince  ourselves  that  all 

G  must  be  equal  to  zero  with  the  exception  of  0  The  con- 
n 

elusion  is  that  the  appropriate  expansion  for  E^,  H  is  the 
following  one : 


E  (a) 
x 


E0xcos<f>0 


sina  -  sin4>0 


+  E. 


V°> 

V*-> 


(3.  12) 


H  (a) 
x 


Hoxcos  <J>0 


+  H 

e 


V  (a) 

e 


(3.  13) 


sin  a-sin<{)0 


%<♦«> 


wherein  E_  ,  Hn  are  the  x-components  of  the  incident  field  and 
x  x 

E  ,  H  two  still  unknown  constants, 
h  e 

In  order  to  evaluate  those  constants,  we  note  that  the 

spectral  components  E  ,  H  as  well  as  E  ,  H  exhibit  simple 

z  z  y  y 

poles  in  the  integrand  (3.  6)  at  values  of  a  given  by 

cos(d+<j>)  =±  - -  (3.  I1*) 

sm80 

When  the  integral  (3-  6)  is  evaluated  in  terms  of  its  singu¬ 
larities,  these  poles  generate  plane  waves  whose  y-dependence 


is  of  the  type: 
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exp(±ky  cosg0)  (3.  15) 

therefore  diverging  for  positive  or  negative  values  of  y. 

In  order  to  have  a  field  solution  which  is  everywhere  boun¬ 
ded,  the  singularities  {3.  1*0  must  be  compensated  forcing 
the  numeratores  in  (3.  10  -  11 )  to  be  zero  at  values  of  a 
given  by  (  3  .  ^*0  •  This  can  be  obtained  by  letting 

E  (a-)  -i C  H  (a-)  =  0  (3.  l6) 

x  p  x  p 

wherein 


cosa  = 
P 


sing0 
and,  consequently: 


sinoO  =±icotgB0 


(3.  IT) 


(3.  18) 


It  can  be  checked  by  ispection  that  relation  (3-  16)  forces 

also  E  ,  H  (see  A.  10  -  11)  to  be  finite  at  , 

y  y  P 

Condition  (3-  l6  )  provides  a  system  of  two  equations  in 

the  two  unknown  E  ,  H  .  Once  these  are  evaluated,  the  spe¬ 
ll  e 

ctral  components  of  the  field  are  completely  known. 


U .  Spectral  field  singularities 


For  a  discussion  about  the  properties  of  the  spectral 


representation  (3.  6)  it  is  convinient  to  deform  the  origi¬ 


nal  integration  contour  T  into  the  two-branch  steepest  de¬ 
scent  path  C  (see  Fig.  5),  so  that 


G(p,<|>)  = 


2iri 


G(  a+(j> )  e  xp  ( -ik  P  sin$0  cosa)da  + 


(1*.  1) 


■E[Residues  of  Gfi^)  in  (T-C)] 


The  function  G(a+<J>)  exhibits  singularities  coincident  with 
those  (3-  9).  Poles  of  the  first  factor  in  (3.  9)  do  accour 
at : 

a  =  -<J>  +  (  — )  n4>  +mr ;  n  =  0  ;  ±1; . (4.  2) 

n 

Poles  of  the  second  factor  do  accour  at 

a+=±(  +  0~ )  +  2nrr  (4.  3) 

If  those  poles  are  located  inside  (T-C),  the  corresponding 
residues  generate  plane-wave  fields  . 

In  particular,  the  pole  o0  corresponds  to  the  incident 
wave;  the  two  poles  a+1  to  the  (geometrically)  reflected  wa¬ 
ves  upon  the  two  half-planes  <p  =  ±n/2  respecti  vely ;  and  the 
poles  a+to  surface  waves  upon  the  same  half-planes. 

Furthermore,  the  integrand  in  (4.  l)  exhibits  saddle 
points  at  a  =  —  4> ± tt  .  In  the  asymptotic  evaluation  of  the 
field,  these  saddle  points  generate  cylindrical  waves  from 
the  discontinuity  line  p  =  0 
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5 .  Surface  vave  incidence 


As  stated  in  Sect.  1,  an  important  problem  is  the  computation 
of  the  surface  wave  conversion  at  the  surface  impedance  discontinuity 
for  a  mixed-path  plane  Earth  propagation.  Accordingly,  we  are  taking 


as  incident  wave  a  TM  (with  respect  to  x)  surface  wave  (see  Fig.  2)  of 
field  x- components  (see  1.  10  and  11 ): 


E1  =  Ec  exp(ikz  cosg0)exp  [-ik  p  sin0o  cos(<)>-4>0)] 

X  X 


H1  =  0 
x 


(5.  l) 


wherein,  to  the  first  order  in  (1  , 


IT 


00=8  ;  <f>0=  --0 


(5.  2) 


The  spectral  x-components  of  the  total  field  are  readily  obtained  from 


(3.  12  and  13),  hence: 

Ec  cos  <p0  Y  (a) 

Eja)  =  — 25 -  +  Eh  — - 

sina  -  sin<j>0  ¥  (<J>0) 


(5.  3) 


H  (a)  *  H 
x  e 


Y  (a) 

e _ 

MO 


(5.  M 


The  fields  for  <£=  -tt/2,  i.e.,  on  the  Earth's  surface  y  <  0  are  of 


interest.  The  three  poles  (see  1+.  2  and  3) 


a  =  -it+0 


a  =  — rr-6 


a-  =  -iT-0 


(5.  5) 
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are  located  near  the  saddle  point  a=  -n  as  far  as  E  (a)  in  concerned, 

x 

while  H  (a)  is  regular  in  the  same  neighbour.  The  two  poles 
x 


a0  =  TT-e 


=  tt+  0 


(5.  6) 


are  located  near  by  the  saddle  point  ci=tt,  again  as  far  as  E  (a)  is  con- 

x 

cerned,  while  H  (a)  is  again  regular  in  the  same  neighbour.  Since  the 
x 

surface  wave  transmitted  in  y  < 0  is  essentialy  the  contribution  of  the 
pole  a  ,  it  follows  that  this  wave  is  still  TM  polarized,  as  it  would 
be  expected,  since  H  (a)  provides  no  contribution  to  it.  For  the  com¬ 
putation  of  this  transmitted  surface  field,  first  we  must  evaluate  the 

constant  E  in  (5.  3).  This  can  be  readily  accomplished  by  forcing  the 
h 


condition  (3.  16),  hence: 

2 

Eov  costp o  Sln  So 

E»- 

1  -  sin  S0  cos  <j>0 


sin  <f>0  +  i  cotgBc 


(5.  T) 


r  (a*)V  (a")  +¥  (qT)'J'  (a+) 

hpep  hpep 

T.  (a+)  T  (a")  -  V(a~)V  (a+) 
hp  ep  hpep 


Now,  in  the  asymptotic  evaluation  of  the  integral  in  (U.  l)  for  large 
kp  sinB0,  only  values  of  the  integrand  near  the  saddle  points  a  =±tt 


are  of  interest,  so  that  a  series  expansion  of  E^(cx-  ir/2)  near  those 
two  points  is  appropriate.  In  particular,  we  get  from  (2.  11  and  12) 


‘K-JT 


-  2b  - 


near 


a=  -  tt. 


E  (-tt+t-  *)  = 
x  2 


E.  sin  6 
x 


+  E 


COST  -  COS0 


0++T 


0  -T 


sin 


~  sin 
2  2 


■  ft+  •  i_±I 
sin  °  sm 


(5-  8) 


A  great  care  should  now  be  taken  in  the  asymptotic  evaluation  of  the 
field;  due  to  the  clustering  of  two  poles  nearby  the  saddle  point.  The 
evaluation  can  be  accomplished  either  by  using  the  Bleistein  method [30] 
or  simply  by  extracting  the  singularity  near  the  saddle  point  [31I .  The¬ 
se  two  methods  are  equivalent  to  the  first  order  and  the  latter  is  most 
convinient  to  our  purposes  since  only  the  term  singular  at  is  here 

of  interest.  This  term,  E'(t),  is  given  by: 

x 


E'(t)  =  lim  sin  — ~  E  (x)  = 
x  2  x 

T  & 


E  sin  6 
x  _ 


..2  0  :  .  2  0  v 
2(sin  ~  -  sm 


+  K 

+  h 


0+-0 


sin9 

sin9 


sm 


(5.  9) 


sm 


0  +  x 


Wherein  the  integration  contour  C  in  depicted  in  Fig.  6.  A  detained 


evaluation  of  the  constant  E  in  terms  of  the  system  narameters  9 

h 

0  ,  B  is  given  in  Appendix  B.  If  not  only  ®  is  small,  but  also 


6 

0  cotg 


«  1 


(9.  10) 


The  result  (5.  13)  is  interesting  inasmuch  as  it  shows  that,  at 
large  distances  from  the  surface  discontinuity,  the  surface  field 
is  transmitted  according  to  the  ususal  tranmission  coefficient  of 
the  refraction  theory.  Furthermore,  the  two  components  of  the  pro¬ 
pagation  constant  are  given  by: 

k  =  k  cos  60  =  k  cos  8  cos  8+  (5.  15) 

Z  -  +  ■  2e+  •  2e~ 

k  =  k  sin  30  cos  0  =  k  sin  B0  cos  0+(l  +  - r — “  ) 

y  2 

Eqs .  (5-  15)  show  that  the  change  in  the  phase  front  direction  is 
negligible  when  the  assumption  (1.  3)  is  made.  However,  expression 

(5-  15)»  at  variance  of  that  for  T,  is  valid  also  when  (1.  3)  is 

.  + 
relaxed, since  only  the  angle  9  is  small. 


S .Excitation  of  the  surface  made  by  a  magnetic  line  source  over  the 


Earth 

In  the  preceeding  section  it  was  shown  that,  in  the  case  of 
highly  conducting  surfaces,  propagation  of  the  surface  wave  is 
poorly  affected  by  the  incidence  angle. 

This  result  suggests  to  examine  the  excitation  of  the  surface 
wave  in  the  two-dimensional  case,  i.e.,  when  the  source  is  taken 
equal  to  a  magnetic  line  source  parallel  to  the  plane  x  =  0,  as 
depicted  in  Fig.  3,  with  n0< <  !•  From  the  analysis  of  Sect.  1 
-see  (1.  13  and  16)  -  we  get  the  formal  expression  for  the  z-compo- 
nent  of  the  total  magnetic  field  along  the  surface  x  *  0,  y  ±  0, 
hence  , s  » 

t  i  ,  f  (6'  11 

H  (p,  -  “)  =  - 2 —  -  j  G(a)  exp  [  ik  pQ  cos(a  -  %)  ]  da 

Z  d  i  O 


wherein  S  is  the  integration  contour  represented  in  Fig.  ^  and 
G(a)  is  given  by  (3.  6  and  9)  with  30  =  90°  and  G^  =1,  hence: 


,  i  'i 
,  cos(  “  -  a) 

G(a)  =  - - - - 

2tti  j 


TT 

4/  (r-  — ) 

h  2  exp(-ikp cos$)a£ 


sin(£  -  “  sin(-|  -  a)  -  a) 


(6.  2) 


The  integration  contour  is  represented  in  Fig.  J  and  can  be  changed 
in  the  alternative  contour  C  provided  that  the  residues  of  the  integrand 
inside  F-C  are  taken  into  account. 


From  result  of  Sect.  U,  it  follows  that  the  pole  location  is  that 


depicted  in  the  same  Fig.  7*  Since  a  is  prescribed  to  run  along  the  curve 


S  of  Fig.  U ,  it  follows  that  the  two  poles  £ 0  and  £  are  always  in¬ 
side  the  closed  contour  f-  0,  for n 0 >  0.  The  corresponding  residues 


since,  in  the  asymptotic  evaluation  of  the  integral  (6.  1)  for  kp  D 
large,  only  values  of  a=T|  0  <<1  are  of  interest.  Note  that  we  used 
result  (2.  11  -  12)  in  order  to  obtain  the  final  expression  (6.  3). 

The  integral  (6.  2)  along  the  two  steepest  descent  paths  C  must 
now  be  computed.  If  : 

/2  k  p  sin  1^°  <<;  1  (6.  4) 


the  position  of  the  poles  £0»£^,£;  ,£  ,  is  sufficiently  remote  from 

the  integration  path  to  produce  any  significant  change  in  the  conventio¬ 
nal  asymptotic  evaluation.  From  another  viewpoint,  decomposition  of 
the  integrand  (6.  2)  in  simple  fractions  and  transformation  via 
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(l.  IT)  in  the  canonical  form  (2.  13)  leads  to  terms  that  are  inversely 
propotional  to  (6.  k),  and  therefore  dominated  by  (6.  3).  This  is  not 
true  for  the  pole  £  .  Expanding  the  integrand  (6.  2)  nearby  £  =  -tt  , 
be  have,  for  a  small. 


cos(—  -  a) 


\U  -  ~) 

h  2 


n 


sin(£  -  sin(~  -  a)  V  2  '  a) 


sm  a 


0  +T  .  0  -T 

sin  -  sm — “ — 

2  2 


0  +Ct  0  +  T 

cos  T  -  cos  a  sin  —  sin  — - - 

2  2 


,  i  =  -TT+T  (6.  5) 


Extracting  the  singularity  at  T  =  -0  and  evaluating  the  integral  as 
in  Sect.  5,  we  get  the  term 


1 

2 


sina  sin  0  sin 


0+-e” 


,.2et  .20,  .  d+  +  a 

(sin  ~  -  sin  -  )  sm  - - 


exp(ikp  cos  0  ) 


(6.  6) 


which  again  dominates  the  remaining  part  of  the  asymptotic  evaluation 


if,  as  we  have  assumed, 
2  - 

kp  sin  0  «  1 


(6.  T) 


The  conclusion  is  that  we  can  take  as  G(at)  the  sum  of  (6.  3)  and 
(6.  6).  For  the  term  (6.  3)  we  have 


1 
*  j 


2  sin  a 


sin  a+  sin  0 


exp(ikp  cos  a  )  exp[  ikpccos(a-n  0}]  da  = 


1 

ir 


r 


2  sin  a  r  ..  .  /  \  , 

-  exp|_  lk  d  cos  (a-  n>  J  da 

sin  a+  sin  0 


S 


(6.  8) 


upon  comparison  with  (1.  lh)  and  when  the  point  P(x,  y )  of  Fig.  3 
lies  on  the  surface  x  *  0. 

By  using  the  same  procedure  as  for  the  integral  (5.  12)  we  get 
for  the  first  term  H  of  the  magentic  field  for  y  <  0 


H  (p, 
Z1 


f>  * 


k  I 


m 


h  Z  \|ir  k  d 


2  7T 

-  exp(-i  r  +  i  k  d)  + 


(6.  9) 


.  - -  .  6  +  r, 

Q-  F  (  '2kd  sin  “ 

+  Ui  sin  ~  exp  [ikd  cos(0  +  r, )} 


f(o) 


When: 


/2k  d  sin 


9  +  n 


»  1 


(  6.  10) 


asymptotic  expansion  of  the  solution  gives 
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H  (p, 
z2  w 


-  *> 


A  1 

m 

2  sin  n 

1 

1 

2 

4  « 

sin  n:'+  sin 

9  \ 

k  d 

exp(-i  ^  +  ikd) 


(6.  11) 


i.e.,  the  cylindrical  wave  radiated  by  the  source  and  its  image  with 
respect  to  the  plane  y  =  0. 

•  .  .  IT 

Similar  results  we  have  for  the  contribution  H  (p,  -  due 

z2  2 

to  the  second  term  (6.  6)  upon  decomposition  in  simple  fraction. 
Computations  are  straightforward. 


7.  Conclusions 

We  presented  a  rigorous  theory  of  propagation  along  a  mixed-path 
plane  Earth.  The  theory  is  based  upon  the  exact  solution  given  by 
Maliuzhinets  for  the  normal  incidence,  properly  extended  to  the  obli¬ 
que  incidence  case.  Maliuzhinets'  solution  is  given  in  terms  of  a 
(relatively)  complicated  function. 

We  have  evaluated  this  function  in  closed  form  for  all  cases  of 
practical  interest. 

First  the  case  of  an  incident  surface  wave  has  been  considered. 
For  highly  conductive  surfaces,  it  was  shown  that  the  transmission 
coefficient  is  real  and  practically  independent  from  the  incidence 
angle.  This  is  not  the  case  when  at  least  one  of  the  surfaces  is 
poorly  conducting.  Although  not  explicitly  examined,  the  solution 


iiMliKHiifi— iiYii  u  •  .v'l 


gives,  in  this  last  case  a  (complex)  angle  depending  trasmission 
coefficient. 

The  direction  of  propagation  of  the  phase  front  is  not  appre¬ 
ciably  changed,  for  the  case  of  highly  conducting  surfaces,  far  from 
the  impedance  discountinuity .  Again,  this  is  not  true  when  at  least 
one  of  the  surfaces  is  poorly  conducting. 

The  case  of  source  excitation  can  be  considered  by  synthesizing 
the  field  from  the  plane  wave  solution.  Examination  of  the  solution 
confirms  previous  results. 


APPENDIX  A 


Relations  between  spectral  components  of  the  field 


For  the  spectral  components  of  the  field  we  have: 


k  x  E  «  oj  u  H 

k  x  H  =  -  w  e  E 

k  •  E  =  0 

k  •  H  =  0 

From  (1.10),  it  follows  that: 


(A.  1) 
(A.  2) 

(A.  3) 
(A.  k) 


k  =  -k  sin  eo  cos  (a+<j>)  i_^  -  k  sin 

=  k  i  +k  i  +k  i 
x  ~x  y  ~y  z  z 

Substituting  (A.  5)  in  (A.  1,  3, 


80  sin(a+(j>) 


i  +  k  cos 

" y 


k)  we  easily  get: 


(A.  5i 


(k2  +  k2)  E  sifv  k  H  -  k  k  E  (A.  6] 

y  z  z  y  x  x  z  x 


Similarly,  by  substituting  (A.  5)  in  (A.  2,  3,  M  we  get: 


2  ~ 

+  k  )  H 
z  z 


-  a.  ek  L  -  k  k  H 
y  x  x  z  x 


(A.  7 


Accordingly: 
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APPENDIX  B 


Evaluation  of  the  constant  E, 


From  (  5.  7 )  the  constant  E  can  be  put  under  the  form: 

h 


E  = 
h 


.  2 

Eox  c°Slto  Sln  S0 


1  -  sin  60  cos  4>e 


sin  $  Q  +  i  cotg  B.  R 


(B.  1) 


¥,  (a  + )  h  (a  )  +  Y  (a + )  /H1  (a  ) 
_  h  p  h  p _ e  p  e  p 

Y  (a+)  /V  (a~)  -  f  (n  +  )  A  (a") 
hp  hp  e  p  ep 


(B.  2) 


wherein,  from  (3-  18) 


±  B 

a  =  i  i  In  cotg 
P  2 


=  l  p 


(B.  3) 


On  the  other  hand,  we  have  from  (2.  U) 


N(a  +x)  =  N(a  )  +  ~  N(a  )  ~ ° 

4TT 


cos  a 


(B.  1*) 


Accordingly,  for  small  values  of  0  , 

2ip)0+  +  (it  -2ip  )0~] 


\(\) 
h  P 


1+-*^  |  (  ti+  2i| 


2  * 


VV  1  +  — ij^  r  ( Tt  -  2ip)9+  +  (ir  +  2ip)e’ 

2  ir  L 


(B.  5) 


Furthermore 


N(z)  =  N(x  +  iy)  = 


lA 

cos  z  exp 


t-j 


exp 


1 

2tt 


fy 


x  +  it 


dt 


cos(x  +  it) 


For  y  »1, 


N(x  +  iy)  N(x  +  i®  ) 
Accordingly,  for  0  small  and 

[e  =ot  §>] 


«  1 


we  have  from  (2.  2) 


y  (a +) 

.  e_.£_ 


y  (a  ) 
e  p 


In  conclusion,  we  have: 


cos 


R  ~ ' 


ip  sin  eo(6  -  6  ) 


E  a 
h 


2 

E0x  cos  <|>0  sin  80 


sin  $0  + 


—  dt 
t 

(B.  6) 

(B.  T) 

(B.  8) 

(B.  < 

(B. 

cotg  B0 


sin  B0  cos  $ 


sin  B0  In  cotg 
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Fig.  3.  Geometry  of  an  uniform  magnetic  source  line  over  a 
mixed-path  Earth. 
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